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Abstract 

We consider a Schrodinger operator with a magnetic field (and no electric 
field) on a domain in the Euclidean space with a compact boundary. We give 
sufficient conditions on the behaviour of the magnetic field near the boundary 
which guarantees essential self-adjointness of this operator. From the physical 
point of view, it means that the quantum particle is confined in the domain by 
the magnetic field. We construct examples in the case where the boundary is 
smooth as well as for polytopes; these examples are highly simplified models 
of what is done for nuclear fusion in tokamacs. We also present some open 
problems. 
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1 Introduction 

1.1 The problem 

Let us consider a particle in a domain Q in R d (d > 2) in the presence of a magnetic 
field B. We will always assume that the topological boundary dVL := f2 \ 51 of SI 
is compact. At the classical level, if the strength of the field tends to infinity as x 
approaches the boundary dfl, we expect that the charged particle is confined and 
never visits the boundary: the Hamiltonian dynamics is complete. At the quantum 
level the fact that the particle never feels the boundary amounts to saying that the 
magnetic field completely determines the motion, so there is no need for boundary 
conditions. At the mathematical level, the problem is to find conditions on the 
behaviour of B(x) as x tends to 9f2 which ensure that the magnetic operator Ha 
is essentially self-adjoint (see Section |2~B]) on C^(Vi) (the space of compactly sup- 
ported smooth functions). These conditions will not depend on the gauge A, but 
only on the field B. One could have called such pairs (f2, ^4) "magnetic bottles" , but 
this denomination is already introduced in the important paper [3] for Schrodinger 
operators with magnetic fields in the whole of M. d having compact resolvents. This 
question may be of technological interest in the construction of tokamacs for the nu- 
clear fusion [30j . The ionised plasma which is heated is confined thanks to magnetic 
fields. 



1.2 Previous works 

The same problem, concerning scalar (electric) potentials, has been intensively stud- 
ied. In the many-dimensional case the basic result appears in a paper by B. Simon 
[24] which generalises results of H. Kalf, J. Walter and U.-V. Schminke (see [15] 
for a general review). Concerning the magnetic potential, the first general result is 
by Ikebe and Kato: in [T3], they prove self-adjointness in the case of = R d for 
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any regular enough magnetic potential. This result was then improved in [251 I26| . 
Concerning domains with boundary, we have not seen results in the purely magnetic 
case. A regularity condition on the direction of the magnetic field was introduced 
in the important paper [2] (Corollary 2.10 p. 853) in order to construct "magnetic 
bottles" in M. d . It was used later in many papers like 

In the recent paper |20j . G. Nenciu and I. Nenciu give an optimal condition for 
essential self-adjointncss on the electric potential near the boundary of a bounded 
smooth domain; they use Agmon-type results on exponential decay of cigenfunctions 
combined with multidimensional Hardy inequalities. 

1.3 Rough description of our results 

As we will see, in the case of a magnetic potential the Agmon-type estimates still 
hold, whereas the Hardy inequalities cannot be used because there is no separation 
between kinetic and potential energy. Actually the point is that we need, to apply 
the strategy of [20] , some lower bound on the magnetic quadratic form Ha associated 
with the magnetic potential A. Our main result is as follows: under some continuity 
assumption on the direction of B(x) at the boundary, for any e > and R > 0, 
there exists a constant C t ^ € K such that the quadratic form Ka satisfies the quite 
optimal bound 



Here | J5(a:) | sp is a suitable norm on the space of bi-linear antisymmetric forms on 
M. d , called the spectral norm. This implies that Ha is essentially self-adjoint if 
there exists 77 > so that |-B(a;)| S p > (1 + rf)D(x)~ 2 where D is the distance to the 
boundary of f2. 

We study then examples in the following cases: 

• The domain f2 is a polytope 

• The boundary d£l is smooth and the Euler characteristic x(cT2) vanishes 
(toroidal domain) 

• The boundary <9f2 is smooth and the Euler characteristic x(dty does not 
vanish (non toroidal domain) 

• The domain is 1R 3 \ and the field is a monopole or a dipolc 

• The domain Q is the unit disk: for any e > and d = 2, we construct 
an example of a non essentially self-adjoint operator Ha with |i?(.T)| sp ~ 
{y/3/2 - e)D{x)~ 2 showing that our bound is rather sharp. 

1.4 Open problems 

The following questions seem to be quite interesting: 

• What are the properties of a classical charged particle in a confining magnetic 
box? Are almost all trajectories not hitting the boundary? 

• What is the optimal constant C in the estimates |-B(x)| sp > CD{x)~ 2 of our 
main result [3~2T From our main results and the example in the unit disk given 
in Section we see that the optimal constant lies in the interval [\/3/2, 1]. 
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2 Definitions and background results 



In this section, we will give precise definitions and related notations. We will also 
review some known results with references to the literature. 

2.1 The domain f2 

In what follows, we will keep the following definitions: Q is an open set in the 
Euclidean space M. d (d > 2) with a compact topological boundary <9f2 = fi \ 0, so 
that either Q or R d \ ft is bounded. 

Definition 2.1 We will denote by dn the distance defined on O by the Riemannian 
metric induced by the Euclidean metric: 

dn{x,y)= inf length^) 

where T x ^ y is the set of smooth curves 7 : [0, 1] — > Cl with 7(0) = x, 7(1) = y. 

We will denote by f2 the metric completion of (O, dn) and by = Q\Q the 
metric boundary ofQ. 

We say that Q is regular if 0^ is compact. 

If fl is regular, dil is compact. In fact the identity map of fl extends to a 
continuous map it from £1 onto Q and ^(Qoa) = dQ. (Q, n) is a "desingularisation" 
of f2. If X = dQ. is a compact C 1 sub-manifold or a compact simplicial complex 
embedded in a piecewise C 1 way, is regular. 

If X = U„£n[0, l]e n with e„ a sequence of unit vectors in M 2 converging to eo, 
then R 2 \ X is not regular, even if d£l = X is compact. 




Figure 1: An example where dX is compact while X^ is not compact 

We will use the following regularity property: 

Definition 2.2 Let us assume that Q is regular. A continuous function f : Q — > C 
is regular at the boundary if it extends by continuity to f2. 

The Lebesgue measure will be |<ia;| and we will denote by (u, v) :— J n uv\dx\ the 
L 2 scalar product and by |ju|| the L 2 norm of u. We will denote by C^°(f2) the 
space of complex- valued smooth functions with compact support in Q. 

2.2 The distance to the boundary 
2.2.1 The distance function 

Definition 2.3 Let us denote by dn the extension of dn by continuity to Q. For x € 
fl, let D(x) be the distance to the boundary Qqo, given by D(x) = min^gn^ dn(x,y). 
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Lemma 2.4 The function D is 1-Lipschitz and almost everywhere differentiable in 
f2. At any point x of differentiability of D, we have \dD(x)\ < 1. 

The inequality \D(x) — D(x')\ < d R (x,x') follows from the triangle inequality for 
d,R. The almost everywhere differentiability of Lipschitz functions is the celebrated 
Theorem of Hans Rademacher [21]; see also [19] p. 65 and [13] . 

2.2.2 Adapted charts for smooth boundaries 

Assuming that the boundary is smooth, we can find, for each point xq € <9f2, a 
diffeomorphism from an open neighbourhood U of xo in M. d onto an open neigh- 
bourhood V of in R^ x , satisfying: 

• xi(F(x)) = D(x) 

• The differential F'(xq) of F is an isometry 

• F(unn) = v n {xi > o}. 

We will call such a chart an adapted chart at the point xq. Such a chart is an 
e— quasi- isometry (see the definition in Section [472]) with e as small as one wants by 
choosing U small enough. 

2.3 Antisymmetric forms 

Let us denote by A^R^ the space of real- valued k-linear antisymmetric forms on the 
Euclidean space R d . The space A 1 R d is the dual of M. d , and it is equipped with the 
natural Euclidean norm: | 5TJj=i a 3 dx 3 \ 2 = 5TJj=i a ]- The space A 2 M d is equipped 
with the spectral norm: if B E A 2 M d , there exists an orthonormal basis of M. d so 
that B = bi 2 dxi A dx 2 + 634^X3 A dx^ + • • • + 6 2 J-i 2d w hh d = [d/2] and bi 2 > 634 > 
• • • > 0; the sequence 612,634, • • • is unique: the eigenvalues of the antisymmetric 
endomorphism B of R d associated with B{x) are ±1612, ±£634, ■ • • , Hb 2 d-i id and 
if d is odd. 

Definition 2.5 We define the spectral norm of B by \B\ sp := 5TJj=i ^2j-i,2j- 

\B\ sp is one half of the trace norm of B, hence |B| sp is a norm on A 2 R d . If d = 2, 
|-B| S p = \B\; if d = 3, |-B| sp is the Euclidean norm of the vector field B associated 
with B, defined by i{B)dx A dy A dz = B where l(B)uj is the inner product of the 
vector field B with the differential form to. 

Remark 2.6 |-B| sp is the infimum of the spectrum of the Schrddinger operator with 
constant magnetic field B in R d . 

2.4 Magnetic fields 

Let us give the basic definitions and notations concerning magnetic fields in a do- 
main fl. The magnetic potential is a smooth real 1-form A on C R d , given by 
A = 57Jj=i Qjdxj, and the associated magnetic field is the 2-form B = dA; more ex- 
plicitly, we have B(x) = J2i< 3< k<d b 3 k(x)dx 3 Adx k with b 3k (x) = d 3 a k (x)~d k aj(x) . 
Let us define now the Schrddinger operator with magnetic field B = dA: 

Definition 2.7 The magnetic connection V = (Vj) is the differential operator 
defined by 

d 



The magnetic Schrodingcr operator Ha is defined by 

3=1 

The magnetic Dirichlct integral Ka = (Ha-\.) is given, for u G Co°(f2), by 




Let us note the commutator formula [Vj,Vfc] = —ibjk which will be very im- 
portant. From the previous definitions and the fact that the formal adjoint of Vj 
is — Vj, it is clear that the operator Ha is symmetric on Cq°(Q.). 

Definition 2.8 We will say that B = dA is a confining field in if Ha is essen- 
tially self-adjoint (see Section ] 2. 6]) . 

2.5 The Riemannian context 

2.5.1 "Regular" Riemannian manifolds 

The context of an Euclidean domain is not the most natural one for our problem. 
In particular, the "regularity assumption" of Definition 12.11 can easily be extended 
to the Riemannian context. Let be a smooth Riemannian manifold. We are 

interested in cases where (H, g) is not complete. Let us recall that g induces on 
Q, a distance d g defined by d g (x,y) = inf 7S r x „ length^) where T XtV is the set of 
smooth paths 7 : [0, 1] — > Q so that 7(0) = x, 7(1) = y. We will denote by f2 the 
metric completion of Q and by = f2 \ O the metric boundary. In the case where 
fl C M. d is equipped with the Euclidean Riemannian metric, f2oo is in general not 
equal to the boundary dCl. 

Definition 12. II is now replaced by: 

Definition 2.9 The Riemannian manifold (Q, g) is regular if 

1. is compact 

2. For any e > 0, every xq £ fioo has a neighbourhood U so that so that U n fi is 
e— quasi-isometric (see Definition ^ ■ 4\ ) to an open set ofM. d with an Euclidean 
metric. 

A function f : Q — > C is regular at the boundary if it extends by continuity to tt. 

2.5.2 Magnetic fields on Riemannian manifolds 

The magnetic potential is a smooth real valued 1-form A on Q, the magnetic field is 
the 2- form B = dA. The norm |2?(x)| sp is calculated with respect to the Euclidean 
metric g Xa . The magnetic potential defines a connection V on the trivial line bundle 
fl x C — > fi by V xf = df(X) — iAf. The magnetic Dirichlet integral is = 
In || V/||g|dx| s where the norm of the 1-form V/(x) is calculated with the dual 
Riemannian norm: ||V/||g = J2ij S^^ft/Va^/ and \dx\ g = 6\dx% ■ ■ ■ dxa\ is the 
Riemannian volume. The magnetic Schrodingcr operator is then defined by: 

H A f = -e- 1 J2^d i (<¥ J 'v a ./) • 

ij 
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2.6 Essential self-adjointness 



In this section, we will review what is an essentially self-adjoint operator and give 
some easy propositions which we were not able to point in the literature. 

2.6.1 Essentially self-adjoint operators 

Let us recall the following 

Definition 2.10 A differential operator P : C^°(f2) — * C^°(fl) is essentially self- 
adjoint in L 2 (f2, \dx\) if P is formally symmetric (for any u,v £ C^°(fl), (Pu\v) = 
(u\Pv)) and the closure of P is self-adjoint. 

A basic criterion for essential self-adjointness is the following (see criterion (4) 
of Theorem X.l and Corollaries in [55] ): 

Proposition 2.11 Let P be as before and formally symmetric. Let us assume 
either that 

(1) there exists E eR so that any solution v G L 2 (Q) of (P — E)v = (in the 
weak sense of Schwartz distributions ) vanishes, 

or that 

(2) there exists X± G C with ±3(A±) > so that any solution v G L 2 (fl) of 
(P — X±)v = (in the weak sense of Schwartz distributions) vanishes. 

Then P is essentially self-adjoint. 

2.6.2 Essential self-adjointness depends only on the boundary behaviour 

Proposition 2.12 Let X be a smooth manifold with a smooth density \dx\. Let 
Lj, j = 1,2 be two formally symmetric elliptic differential operators of degree m on 
L 2 (X, \dx\) and let us assume that L\ is essentially self-adjoint and Li — L\ = AI 
is compactly supported. Then L2 is essentially self-adjoint. 

Proof. - It is enough to show that L2 — ci is invertible for c real and large enough. We 
have L2 — ci = (id + AI(Li — ci)^ 1 ) (L\ — ci). Moreover the domain of L\ contains 
.ff™ (the space of compactly supported H m functions). So that ||A/(ii — c« ) 1 1 j = 
Ofc- 1 ). □ 

2.6.3 Essential self-adjointness is independent of the choice of a gauge 

Proposition 2.13 Let X be a smooth manifold with a smooth density \dx\. Let 
us consider a Schrddinger operator and A2 = A\ + dF with F G C°°(X, R). 
Then, if Ha 1 is essentially self-adjoint, Ha 2 is also essentially self-adjoint. 

Proof - We have formally (as differential operators) 

H A2 =e lF H Al e~* F . 

Hence, Ha 2 — ci = e lF {LLax —ci)e~ lF . The domain D2 of the closure of Ha 2 
(defined on C§°(X)) is e lF times the domain D\ of the closure of Ha x - The re- 
sult follows from the fact that e ±lF is invertible in L 2 and an isomorphism of the 
domains. □ 



3 Main results 

Let us take Ha with domain T>(Ha) = C(f(Q). As explained in the introduction, 
we are looking for growth assumptions on |-B| sp close to dVL ensuring essential self- 
adjointness of Ha- We formulate now our main results: 
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Theorem 3.1 Let us take d — 2. Assume that dQ is compact with a finite number 
of connected components and that B(x) satisfies near dfl 

\B(x)\ sp >(D(x))- 2 , (3.1) 

then the Schrodinger operator Ha is essentially self-adjoint. This still holds true 
for any gauge A 1 such that dA' = dA = B. 

Theorem 3.2 Let us take d > 2. Assume that Q is regular and that there exists 
n > such that B(x) satisfies near <9f2 

\B(x)\ sp > (l + n)(D(x))- 2 , (3.2) 

and that the functions 

n ^ x) = \mu P (3 - 3) 

are regular at the boundary f^oo (for any 1 < j ' < k < d) (see Deflnition \2.ty) . Then 
the Schrodinger operator Ha is essentially self-adjoint. This still holds true for any 
gauge A' such that dA' = dA = B. 

Remark 3.3 If fl is defined (locally or globally) by £1 := {x E R d | f(x) > 0} with 
f : R d -> E smooth, df{y) ^ for y G dtt, then f {x) - \df{x)\D{x) for x close to 
dQ. And we can replace in the estimates \3.'0j) D(x) by f(x)/\df(x)\. 

Theorem 13 . 2 1 can be extended to Riemannian manifolds as follows: 

Theorem 3.4 Let (f2,g) be a regular Riemannian manifold with a magnetic field 
B = dA. Let us assume that ||-B|| sp > (1 + e)D~ 2 near f2oo and that, for each 
xq G tloo, the direction n{x) of B, calculated with the metric g XQ (i.e. using the 
trivialisation of the tangent bundle associated with g Xo ), has a limit as x — > xq, then 
Ha is essentially self-adjoint on Cg°(f2). 

The exponent 2 of the leading term in Equations (|3.ip and (|3.2[) is optimal, as 
shown in the following 

Proposition 3.5 For any < a < v3/2, there exists a magnetic field B such that 
Ha (with dA = B) is not essentially self-adjoint and such that the growth of \B\ sp 
near the boundary dfl satisfies 

We prove this proposition in Section [STU in the case d = 2, but the proof can be 
easily generalised to larger dimensions. 

As a consequence of this proposition, together with Theorem 13.11 (respectively 
13.21 ). we get that the optimal constant in front of the leading term (D(x))~ 2 is in 
h/3/2,1]. 

Hence we see that the situation for confining magnetic fields is not the same 
as for confining potentials (for which the optimal constant is 3/4, hence is smaller 
than V3/2). 

Indeed this is due to the difference between the Hardy inequalities in the two 
situations: the term 1/(4D 2 ) does not appear in the magnetic case, as it does in 
the case of a scalar potential, where it plays the role of an " additional barrier" . 

4 Proof of the main results 

In this Section, we prove Theorems 13. li [3~2l and [34l using the method of [20] which 
we first review. 
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4.1 Agmon estimates 



The following statement is proved, using Agmon estimates pQ, in |20| : 

Theorem 4.1 Assume that d£l is compact, and that there exists c £ R such that, 
for all u e C§°(fl), h^iu) — J Q D(x)~ 2 \u(x)\ 2 \dx\ > c||uj| 2 . Then Ha is essentially 
self-adjoint. 

Reading the proof in [2D], one sees that the only property of fl which is used is 
that the function D(x) is smooth near the boundary and satisfies |dD(x)| < 1- One 
can extend the proof to the case where dfl is not a smooth manifold by using the 
properties of the function D described in Lemma 12.41 The fact that is bounded 
does not play an important role, only the compactness of dQ is important. The 
essential self-adjointness of Ha results from the Proposition ^. Ill and the following 

Theorem 4.2 Let v G L 2 (fl) be a weak solution of (Ha — E)v = 0. Let us assume 
that there exists a constant c > such that, for all u £ Cq°(0), 

(u\(H A - E)u) f JTT^\ dx \ 9 ^ C IH! 2 • 

J{xen | D(x)<i} u \ x ) 

Then v = 0. 

For the reader's convenience, we give here the proof of Theorem 14.21 following the 
strategy of [20] in a slightly simplified way. 

Proof. - The proof is based on the following simple identity ([20]) 

Lemma 4.3 Let v be a weak solution of (Ha — E)v = 0, and let f be a real-valued 
Lipschitz function with compact support. Then 

(fv\(H A -E)(fv)) = (v | \df(x)\ 2 v) . (4.2) 

Let us give two numbers p and R satisfying respectively < p < | and 1 < R < 
+oo. We will apply identity (|4.2p with / = F(D) where F(u) the piecewise smooth 
function defined by 

for u< p and for u > R + 1 
2(u — p) for p < u < 2p 

F(u) = { u for 1p < u < 1 

1 for 1 < u < R 
R + 1 - u for R < u < R + 1 




Figure 2: The function F 

We have \df\ 2 = F'(D) 2 almost everywhere. From the inequality (|4.1[) applied 
to fv, we get: 

((H A E)(fv) \fv)> ( \v\ 2 \dx\ g + c\\fv\\ 2 . (4.3) 

J2p<D(x)<l 
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On the other hand, using the explicit values of df and Equation (|4.2[) . we get: 

((H A -E)(fv)\fv)< ^ p < D(x) < 2p \v?\dx\ a + --- 
■ ■ ■ kp<D{x)<M 2 \ dx \s+ Ir<d{x)<r+i \ v \ 2 \ dx \g ■ 

Putting together the inequalities (|4.3[) and (|4.4[) . we get 

c\\fv\\ 2 <4[ \v\ 2 \dx\ g + [ \v\ 2 \dx\ g . (4.5) 

J p<D(x)<2p JR<D{x)<R+l 

Taking p — > and i? — > +00 in the inequalities (|4.5p , we get that the L 2 norm of f 
vanishes. 

□ 

4.2 Quasi-isometries 

In section 5 we give examples which have smooth boundaries (excepting the convex 
polyhcdra fsection l5.ip ). In order to build new examples, like non convex polyhcdra, 
one can use quasi-isometries. 

Definition 4.4 Given < c < C, a (c,C)-quasi-isometry of VL\ onto £l 2 is an 
homeomorphism of F : Qi onto il 2 whose restriction to fli is a smooth diffeomor- 
phism onto Q 2 and such that 

Vx,y € f2i, cd R (x,y) < d R (F(x), F(y)) < Cd R (x,y) . 

An e— quasi- is ometry is an (1 — e, 1 + e) quasi-isometry. 

Lemma 4.5 We have the bounds 

\\F'\\ < C, UF-^'W < c-\ \det(F')\ < C d , cD 1 (x) < D 2 (F(x)) < CD^x), 
where, for i = 1,2, Di(x) denotes, for any x <E f2j, the distances to the boundary 

We will start with a magnetic potential A 2 in and define A\ = F*(A 2 ). We want 
to compare the magnetic quadratic forms h A2 (u) and h Al {u o F) as well as the L 2 
norms. We get: 



Theorem 4.6 Assuming that, for any u € Cq°(C1 2 ), 

J! 

we have, for any v <E Cq°(£Ii), 



f \u\ 2 

h A2 {u)>K LL. \dx 2 \-L\\u\ 

JClo u 2 



h Al (v)>K^) d+2 jf ^dx^-Lc 2 ^ 



In other words, we can check that H Al is essentially self-adjoint from an estimate 
for h A2 using Theorem 14. II 

Proof - Let us start making the change of variables x 2 = F(x\) in the integral 
h A2 (u). Putting v = uo F, we get h A2 (u) = J n± Al v(xi)\\ 2 g \det{F' (xi))\\dx!\ 
where g is the inverse of the pull-back of the Euclidean metric by F. Using Lemma 
14.51 we get the estimate. □ 
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4.3 Lower bounds for the magnetic Dirichlet integrals 
4.3.1 Basic magnetic estimates 
Lemma 4.7 For any u G Cq°(0), we have 

h A (u) > \(b 12 u\u)\ + \(b 34 u\u)\ H h \(b 2 d-i,2d u \ u )\ ■ 

Proof. - We have 

\(b 12 u\u)\^\([V ll V 2 ]u\u)\<2\(V 1 u\V 2 u)\< [ (\V lU { 2 + \V 2 u\ 2 )\dx\ 



We take the sum of similar inequalities replacing the indices (1, 2) by (3, 4), • ■ ■ , (2d— 
1,23). □ 

Lemma 4.8 Let fl be a regular open set in M. d . Let xq £ fioo and assume that 
B(x) does not vanish near the point Xq and that the direction of B is regular near 
xq. Let A be a local potential for B near xq, then, for any e > 0, there exists a 
neighbourhood U of Xq in M. d so that, for any (f> E Cfi°(U D £1), 



> (1 - c) / IBMI-pMaOndzl . (4-6) 
Ju 

where \B(x)\ sp is defined in Definition \2.5\ 

Proof. Let us choose U so that, for all x G U CI, \n(x) — n-(x )|Euci < € \J d(d-i) ' 
where | J2i<j a-ijdxi A dxj | EucI = X)i<j a f?> by applying Definition 12. 21 to n(x) at the 
point xq- We choose orthonormal coordinates in R d so that n(xo) = n\ 2 dx\ /\dx 2 + 
ris&dxs A dx± + • • • with ri2k-i,2k > and Ti2fc— i,2fc = 1- From Lemma 14.71 we 
have, for e C* °°(fi n [/), 

M0)> / \B(x)\ sp (n 12 {x)+n 34 {x) + ---)\<j){x)\ 2 \dx\ 
Ju 

and ni 2 (x)+ri34(x) + - • • > 1 — e, because the Euclidean norm of n(x) is independent 
of the orthonormal basis. □ 

Remark 4.9 The estimate is optimal in view of Remark \2.6l 

4.3.2 The 2-dimensional case 

Theorem 4.10 Let us assume that dfl C B(0, R) and that dtt has a finite number 
of connected components. If d = 2 and if B does not vanish near dfl, then there 
exists c R eR so that, Vu G Cg°(Q,), 

h A (u) > I \B\\u\ 2 \dx\ - c R \\u\\ 2 . (4.7) 

JunB(0,R) 

Proof- As B does not vanish near dCl, the sign of B is constant near each connected 
component of dfl. Let us write fl C U^ =1 0; with f2; open sets such that fliDdfl — 0, 
B > on £l 2 and B < on SI3. We can assume that fl 2 and O3 are bounded. Take a 
partition of unity j = 1, 2, 3, so that, for j = 2, 3, <frj G C{f(Slj), and 4> 2 j = 1- 
Now we use the IMS formula (see [23] ) 

h A (u) = Y, h A(M) - f fel#*N M' \dx\ • (4.8) 

(=0 ^° \l=0 / 



with the lower bound of Lemma 14771 in fli n Q for I = 2, 3 and the lower bound for 
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4.3.3 The case d > 2 



Theorem 4.11 Let us assume that dtt C B(0, R). Assume that B = dA does not 
vanish near dtt and that the functions n,jk(x) are regular at the boundary dfl, then, 
for any e > 0, there exists C e ^n > so that, Vit G Cq°(Q), 

h A (u) > (1 - e) / \B\ sp \u\ 2 \dx\ ~ C t . R [ \u\ 2 \dx\ . (4.9) 

JnnB(o,R) Jn 

Proof - We first choose a finite covering of floo by open sets Ui, I = 1, • • • N of M. d 
which satisfies the estimates of Lemma 14.81 We choose then a partition of unity 
<pl,l = 0,--- ,N with 

• For/>l, <f>i€CF(Ui) 

• 0o is c^(n) 

• Ez <\>\ = 1 in n 

. sup£,|#;| 2 = M. 

Using the estimates given in Lemma T4. 81 for I > 1 and the fact that J^i \d<t>i\ 2 is 
bounded by M, we get, using IMS identity (|4.8p . the inequality ([S^]) . 

□ 

4.4 End of the proof of the main theorems 

Using Theorem 14. 1[ it is enough to show that there exists c G M such that, for all 
u G C °°(tt), 

&a(«)> / |Z?(a;)r 2 | U (x)| 2 |dx|- C || U || 2 , 

Jnr\B(0,R) 

under the assumptions of Theorems 13.11 and 13.21 This is a consequence of Theorem 
OO] for d = 2 and Theorem HH] for d > 2. 

The proof of Theorem [33] is an adaptation of the case of an Euclidean domain. 
The partition of unity is constructed using only the distance function which has 
enough regularity. We use also the fact that near each point xq of the boundary 
the metric is quasi- isometricaHy close to the Euclidean metric g Xo . 

5 Examples 
5.1 Polytopes 

A polytope is a convex compact polyhedron. Let f2 be a polytope given by 

n = ng^x | Li(x) < o} , 

where the Lj's are the affine real- valued functions 

d 

Li(x) = n ijXj + Oj . 

3=1 

We will assume that, for i = 1, • • • , d, Ej=i n lj = 1 (normalisation) and n,i ^ 
(this last condition can always be satisfied by moving Q by a generic isometry) . We 
have the 
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Theorem 5.1 The operator Ha in f2 with 



A 



tjiiLi n 2 iL 2 



1 1 



+ 



) 



dX2 , 



is essentially self-adjoint. 



Proof. - We have 




+ • • • J dx\ A dx2 + bjdxj A dx 



2 • 



and Z) = mini<i<jv \Li\. So that Z? = b^dxi A + j=z^ &dx j A with 
&12 > -D~ 2 . We then apply directly Lemma l4~7l and Theorem 14. 11 



5.2 Examples in domains whose Euler characteristic of the 
boundary vanishes ("toroidal domains"). 

Let us assume that dQ is a smooth compact manifold of co-dimension 1 and denote 
by j : dtt — ► R d the injection of 9f2 into K d . A famous theorem of H. Hopf (see 
[5J [T2]) asserts that there exists a nowhere vanishing tangent vector field to <9£1 (or 
1-form) if and only if the Euler characteristic of dfl vanishes. 

Theorem 5.2 Let us assume that the Euler characteristic of d£l vanishes (we say 
that f2 is toroidal). Let Aq be a smooth Inform on f2 so that the 1—form on dfl 
defined by uj = j*(Ao) does not vanish, and consider a Inform A in f2 defined, near 
dfl, by A — Ao/D a . We assume that either a > 1, or a = 1 with the additional 
condition that for any y £ dfl, > 1. Then Ha is essentially self-adjoint. 

Remark 5.3 The existence of oj is provided by the topological assumption on dQ. 
This works if f2 C K 3 is bounded by a 2-torus. It is the case for tokamacs. 

Proof. - We will apply Theorem 13.21 We have to check: 

• The uniform continuity of the direction of the magnetic field or the extension 
by continuity to It has to be checked locally near the boundary DO.. We 
will use an adapted chart (see section [2. 2. 2[) . 

In these local coordinates we write Aq = a\dxi+{3 with [3 = a 2 dx2 + - ■ -+addxd 
and <jj = 02(0, x')dx 2 + • • • + a<j(0, x')dxd so we get 



Thus we get that the direction of B is equivalent as x\ — > + to that of dxiAu 
which is non vanishing and continuous on Q. 

• The lower bound (|3.2[) \B\ sp > (1 + if)D~ 2 near dtt. The norm of B near the 
boundary is given, as x — > y by 



□ 




\B{x)\ 



a\u{y)\/D 



,a+l 



Therefore we conclude that the hypotheses of Theorem 13.21 are fulfilled. 



□ 
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Remark 5.4 The asymptotic behaviour of B(x) as x — > dQ is 

adx! Aw(y) 

w f^+l — • 

If follows that lo and a depend only of B and are invariant by any gauge transform 
in n. 

Remark 5.5 If d = 3, the magnetic field B can be identified with a vector field B 
in n defined by 

i (^ifj dxi A dx2 A dx% = B 

as in Section \2.3\ Using the induced Riemannian structure, we can identify any 
1-form lo on 90 with a vector field Co. Moreover 90 is oriented by any 2-form 
£ = b(u)dx\ A dx2 A dx^ with v any outgoing vector field near 90. Using the 
previous identifications, the asymptotic behaviour of B near 90 is given by 

B ~ ar (lo) /D a+1 , 

where r is the rotation by +ir/2 in the tangent space to 90. 

It means that B is very large near 90 and parallel to 90. From the point of view 
of classical mechanics, the trajectories of the charged particle are spiralling around 
the field lines and do not cross the boundary. It would be nice to have a precise 
statement. 



5.3 Non toroidal domains 
5.3.1 Statement of results 

We try to follow the same strategy than in Section 15.21 but now we will allow the 
1-form lo on X = 90 to have some zeroes. This is forced by the topology if the 
Eulcr characteristic of 90 does not vanish. We need the 

Definition 5.6 A 1-form to on a compact manifold X is generic if to has a finite 
number of zeroes and dio does not vanish at the zeroes of lo . 

We have the 

Theorem 5.7 Let O C W 1 with a smooth compact boundary X = 90. Let Aq be a 
smooth 1-form in M. d so that lo = j^(Ao) is generic. We assume also that, at each 
zero m of lo, 

\dLo(m)\ sp > 1 , (5.1) 

where the norm |dw(m)| sp is calculated in the space of anti- symmetric bi-linear 
forms on the tangent space T m dVL. Then, if A is a 1-form in O such that near X, 
A = Aq/D 2 , B = dA is confining in O. 

We see that the field need to be more singular than in the toroidal case. We 
could have taken this highly singular part only near the zeroes of lo. 



5.3.2 Local model 

We will work in an adapted chart at a zero of lo. We take A = Aq/xi 2 with 
j*(Ao) — lo, we have: A = a\dx\ + [3 and (3(0) = 0. 
We have 

B = —~ + dx\ A p + 0(x7 1 ). 
x{ 



14 



Applying the basic estimates of Lcmma l4.7l in some orthonormal coordinates in R d_1 
so that <iw(0) = 623^2 A dx% + • ■ • , we see, using the assumption (|5.1|) . that there 
exists a neighbourhood U of the origin and an 77 > so that, for any u £ C^°(/7), 

h A (u)>(l + v) I Hr\dx\ . 



5.3.3 Globalisation 

Near each zero of w, we take a local chart of M. d where A is given by the local model. 
Such a chart is an e-quasi-isometry (see 14.41 ) with e as small as one wants. This 
gives the local estimate near the zeroes of u). The local estimate outside the zeroes 
of oj is clear because we have then |-B| sp > C/D 3 with C > 0: this follows from the 
estimates in Section I5~2l with a = 2. We finish the proof of Theorem 15.71 with IMS 
formula and the local estimates needed in Theorem 14. II 



5.4 An example of a non essentially self-adjoint Schrodinger 
operator with large magnetic field near the boundary 

Let us consider the 1-form defined on f2 = {(x,y) £ R 2 | x 2 + y 2 = r 2 < 1} by 
A = a(xdy — ydx) / (r — 1) where < a < y/3/2. The magnetic potential A is 
invariant by rotations. Then 

Theorem 5.8 The operator Ha is not essentially self-adjoint. 

The corresponding magnetic field B writes B(x,y) = ^fS-^l dx A dy , and, near 
the boundary, ~ a/(D(x)) 2 . We have, in polar coordinates (r, 9), 

^ d 2 1 d 2iar d a 2 r 2 



dr 2 r dr r — 1 88 (r — l) 2 

Hence the operator Ha splits as a sum Ylmez, ^A,m where Ha,™ acts on functions 
e' f(r). We will look at the m = component: Theorem 15.81 follows from the 

Lemma 5.9 If < a < V3/2, on the Hilbert space L 2 (]0,l[,rdr), the operator 

_ d 2 Id a 2 r 2 
dr 2 r dr (r — l) 2 

is in the limit circle case near r = 1 and hence is not essentially self-adjoint. 

Proof.- Let U be the unitary transform U : u — > r x / 2 u from L 2 (]0, 1[, rdr) onto 
£ 2 (]0, 1[, dr). Then K = UHIJ- 1 is given by 

d 2 1 a 2 r 2 



dr 2 4r 2 (r - l) 2 
K is known to be in the limit circle case at r = 1 (Theorem X.10 in [25].) □ 

5.5 Singular points 
5.5.1 Monopoles 

We will first discuss the case of monopoles in K 3 . Here is R 3 \ 0. 
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Definition 5.10 The monopolc of degree m, m G Z \ 0, is the magnetic field 
B m = (m/2)p*(er) where p : K 3 \0 — > § 2 is the radial projection and a the area form 
on 8 2 . In coordinates 

to xdy Adz + ydz A dx + zdx A dy 

2 (3.2 + y 2 + z 2)3/2 

Remark 5.11 Lei us noie, /or comparisons with the case where dfl is of codimen- 
sion 1, that \B m \ sp = ^p-r~ 2 . 

The flux of B m through S 2 is equal to 27tto. This is a well-known quantisation 
condition which is needed in order to build a quantum monopole. In order to define 
the Schrodinger operator H m , we first introduce an Hcrmitian complex line bundle 
L m with an Hcrmitian connexion V m on £1 with curvature B m . We first construct 
L m and V m on § 2 and then take their pull-backs: V m in a direction tangent to a 
sphere is the same and V m vanishes on radial directions. We have, using spherical 
coordinates, 

_d^_2d_ 1 „ 
c/r^ r or r z 

where K m is the angular Schrodinger operator on S 2 (discussed for example in 
]27|). Let us denote by A™ the lowest eigenvalue of K m . The self-adjointness of H m 
depends of the value of A™. As a consequence of Weyl's theory for Sturm-Liouville 
equations, H m is essentially self-adjoint if and only if A™ > 3/4. From p"7l [T8l [27] 
(sketched in Section f5. 5. 2| . we know that A™ = \m\/2 so that 

Theorem 5.12 The Schrodinger operator H m (monopole of degree m) is essen- 
tially self-adjoint if and only if \m\ > 2. 



5.5.2 The spectra of the operators K m , the "spherical Landau levels" 

These spectra are computed in [TTJ [T5] and in the PhD thesis [57] . We sketch here 
the calculus. Recall that K m is the Schrodinger operator with magnetic field ma/2 
where a is the area form on S 2 . The metric is the usual Riemannian metric on S 2 : 

Theorem 5.13 The spectrum of K m is the sequence 

A fe = - (k(k + 2) - m 2 ) , k = \m\, \m\ + 2, • ■ ■ , 

with multiplicities k+1. In particular, the ground state A| m | of K m is \m\/2, with 
multiplicity \m\ + 1. The ground state is exactly the norm of the magnetic field. 

If m = 0, the reader will recognise the spectrum of the Laplace operator on § 2 . 

We start with the sphere § 3 with the canonical metric. Looking at S 3 C C 2 , 
we get an free isometric action of Sg on S 3 : 9. (2:1,2:2) = e l9 (21,22)- The quotient 
manifold is § 2 with 1/4 times the canonical metric; the volume 27r 2 of S 3 divided 
by 2n is n which is one forth of 4tt. 

The quotient map S 3 — ► S 2 is the Hopf fibration, a S 1 — principal bundle. The 
sections of L m over S 2 are identified with the functions on § 3 which satisfy f{9z) = 
e lm9 f{z). With this identification of the sections of L m , we have 

K m = - (A S 3 - to 2 ) , 

where 1/4 comes from the fact that the quotient metric is 1/4 of the canonical one 
and to 2 from the action of d@ which has to be removed. It is enough then to look 
at the spectral decomposition of Ag3 using spherical harmonics: the kth eigenspace 
of A§3 is of dimension (k + 1) 2 and splits into k + 1 subspaces of dimension k+1 
corresponding to m = —k, —k + 2, ■ • ■ , k. 



1G 



5.5.3 A general result for fi = R d \ 

In this section = M d \ and B is singular at the origin. 

Theorem 5.14 // lim^^o M^-^MIsp = +°° and, for any x ^ 0, the direction 
n(tx) has a limit as t — > + , then Mb is essentially self-adjoint 



Proof. - The proof is essentially the same as the proof of Theorem 13 . 21 except that in 
the application of IMS method, we have to take a conical partition of unity whose 
gradients can only be bounded by □ 

5.5.4 Multipoles 

Let us denote, for x G R 3 , B x the monopole with centre x: B x = t*(£?2) with t x 
the translation by x and B^ the monopole with m = 2. If P (^) is a homogeneous 
linear differential operator of degree n on M 3 with constant coefficients, we define 
Bp = P(B x ) x= o- Then Bp is called a multipole of degree n. All multipoles are 
exact! It is a consequence of the famous Cartan's formula: if P is of degree 1, hence 
a constant vector field, 

B v = £vB = d(i(V)B ) . 

A multipole of degree 1 is called a dipole; viewed from very far away, the magnetic 
field of the earth looks like a dipolc. 

Theorem 5.15 If By = dAy is a dipole , Ha v is essentially self-adjoint. 

Proof. - Because By is homogeneous of degree —a = —3, it is enough, using 1 5.141 to 
show that By does not vanish. V is a constant vector field, hence up to a dilatation, 
we can take V = d/dz. We have 

d xdy Adz + ydz hdx + (z — t)dx A dy 

9/92 = di\t=0 {x 2 +y 2 + (z _ t) 2f2 ' 

which gives 

_ 3xzdy Adz + Zyzdz A dx + (2z 2 - x 2 - y 2 )dx A dy 
9/dz = (x 2 +y 2 + z 2 f/ 2 " 

The form Bq/q z does not vanish in f2. □ 

Remark 5.16 We do not know if all multipoles of degree > 2 are essentially self- 
adjoint. 
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